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ABSTRACT 

Spheroidal components of spiral galaxies have been considered the only 
dynamically important component in gravitational lensing studies thus far. Here 
we point out that including the disk component can have a signihcant effect, 
depending on the disk inclination, on a variety of lensing properties that are 
relevant to present studies and future surveys. As an example, we look at the 
multiple image system B1600+434, recently identihed as being lensed by a spiral 
galaxy. We hnd that including the disk component one can understand the 
fairly large image separation as being due to the inclination of a typical spiral, 
rather than the presence of a very massive halo. The fairly low magnihcation 
ratio can also be readily understood if the disk is included. We also discuss how 
such lensed systems might allow one to constrain parameters of spiral galaxies 
such as a disk-to-halo mass ratio, and disk mass scale length. Another example 
we consider is the quasar multiple-lensing cross section, which we hnd can 
increase many-fold at high inclination for a typical spiral. Finally, we discuss 
the changes in the gravitational lensing effects on damped Lyman alpha systems 
(DLAS) when disk lensing is included. 


Subject headings: gravitational lensing—galaxies:structure—galaxies:spiral 


quasars:absorption lines 
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1. Introduction 

Galaxies responsible for gravitational lensing have typically been modeled as one- 
component systems in which the entire mass distribution of the galaxy lens, visible plus dark 
matter, is fairly well approximated by a singular isothermal sphere (SIS): p{r) = Po{ro/rY 
(see e.g. Narayan and Bartelmann 1997). This is probably justihed for E/SO galaxies where 
both components have a surface density fairly independent of orientation (a small ellipticity 
is required only to understand the presence of four-image systems) and the high central 
density represents the high density of the visible matter. By contrast, in spiral galaxies the 
disk component has a very different structure from the dark matter halo and the bulge, 
which results in signihcantly different projected surface density depending on orientation. 
Given that the projected surface density of the disk changes by a factor of radius/thickness 
~ 15 between a face-on and an edge-on spiral, and that the dark and visible components are 
known to contribute roughly equally to the total mass within 4-5 disk scale lengths, one 
would expect signihcant effects on lensing properties for images within a couple of optical 
radii. 

The identihcation of the lensing galaxy in the double image system B1600-I-434 as 
a nearly edge on spiral by Jaunsen and Hjorth (1997) gives us a perfect example to 
demonstrate the importance of including the disk component in modeling lensing by spirals. 
However, an understanding of the lensing geometry of spirals is important for several 
reasons besides that of modeling observed systems. First, recently the gravitational lensing 
effects on damped Lyman alpha systems (DLAS) have been investigated (Smette, Glaeskens, 
and Surdej 1997; Bartelmann and Loeb 1996). At least some DLAS are believed to be due 
to absorption in protogalactic disks or spiral galaxies (Wolfe 1988, 1997), and if the gas 
column density falls off exponentially, one is necessarily dealing with lensing within a few 
exponential scale lengths of the galaxy center, which is where we expect disk lensing to be 



4 


important. Secondly, the extra parameters added by introducing a disk component are the 
physically interesting ones of disk surface density and scale length, so that the observational 
study of gravitational lensing by spirals will provide new constraints on the distribution of 
mass in the various components of spirals. Although only a few such systems are currently 
known, surveys such as the Sloan Digital Sky Survey will catalog a large number of quasars, 
among which there should be a signihcant number lensed by spiral galaxies. This is why it 
is also important to understand the effect of disk lensing on statistical properties such as 
the cross section for multiple-image lensing by spirals. 

In this paper we hrst explain our method of modeling spiral galaxy lenses as singular 
isothermal spheres plus inhnitesimally thin disks (SIS-|-disk). We start with the simple 
example of a constant density disk, to show the effect of adding a disk. We show that such 
a model fits many of the observations on B1600-I-434, and discuss the results of modeling, 
more properly, with an exponential disk. We then discuss the dependence of the lens cross 
section for multiple lensing of background sources on the inclination of the disk, and discuss 
also how it would enhance the gravitational lensing effects on DLAS. We hnish with the 
conclusions and implications from our work. 


2. The Model 

We model a spiral galaxy as consisting of two mass components: a spherical halo and 
an inhnitesimally thin disk. We assume that the halo has the density prohle of a singular 
isothermal sphere, while the disk is treated as having surface density that is either constant, 
or falls off exponentially. For simplicity and clarity we start by treating the case of the 
constant density disk. 


The orientation of the disk dehnes the coordinate system, such that the x and y 
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directions are along the major and minor axis, respectively, of the projected disk, and the 
disk center defines the origin. The model has then 6 parameters: the position of the sonrce 
{Px,Py), the inclination of the disk 70 , the velocity dispersion of the SIS (the SIS can 
be eqnivalently be described by p(r) = a'^/2Gnr'^), and the central snrface density So and 
characteristic length of the disk. For a constant density disk this characteristic length is 
the radins of the disk, while in the exponential disk case this becomes the scale length, 
Rg. Thus, measuring the position of the lens, the inclination of the disk, the two image 
positions {6x, Oy), and the magnihcation ratio completely determines the system. 


The inclined circular disk projects onto an ellipse. We denote the axis ratio of an 
ellipse in general as 7 = h/a, and that of our disk as 70 . Any point {9x, Oy) external to the 
ellipse will dehne a new confocal ellipse with axis ratio 7 '. Then the reduced deflection 
angle, modihed from Schramm (1990), is given by 


ay 




2Er 


1 -y 

l-7o' 


2So, 1-7' 

Sc. ^'(l-7o^)‘ 


( 1 ) 

( 2 ) 


Here is the usual critical density dehned as 

AttGDiDis 


= 




( 3 ) 


and Di is the lens, source, or lens-to-source angular diameter distance. If the point {Ox, Oy) 
is inside the disk, then 7 ' is just replaced with 70 . 


3. The System 

In the gravitational lens system B1600+434 identihed by Jackson et al. (1995), the 
lensed quasar is at a redshift Zg = 1.61, the image separation is V'A, and the intensity ratio 
is Ia/Ib = 1-3 in the radio (in the optical, image B suffers significant extinction). The 
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lens galaxy, identified as a nearly edge on spiral by Jaunsen and Hjortli (1997), has an 
estimated photometric redshift zi ~ 0.4, and axis ratio a/h = 2.4 ± 0.2, or 7 o = 0.42. (It 
is possible that this galaxy is actnally an SO, as snggested to us by C. Kochanek, private 
communication, Jannuary 1997, but this will not effect our analysis as there is still a 
disk.) To emphasize that this paper is meant as a qualitative exploration of the SIS+disk 
model, and because zi is still uncertain, we treat all measurements as exact and include no 
estimates of errors. The positions of the images in the lens-centered, axes-dehned coordinate 
system are given in Table I. Note that the lensing galaxy is just 0735 from image B. In all of 
our calculations, we assume hi = 1, go = 0.5 and Ho = 100hkms“^Mpc“^. Also we calculate 
angular diameter distances assuming all the matter in the universe is smoothly distributed 
(see Schneider, Ehlers and Falco 1992). 

Treating the lens as a singular isothermal sphere requires a velocity dispersion of 200 
km/s to account for the image separation. Of course, in this case the lens must be collinear 
with the images, which is not the case observationally (one would need to introduce some 
ellipticity in the dark matter distribution to avoid this). However, we can obtain an estimate 
of the implied relative magnihcation and the time delay by shifting the lens to the collinear 
position along a line perpendicular to AB. The implied time delay is AtAB = 24/i“^ days, 
and the magnihcation ratio is 3.8. (The results are similar if the lens position is shifted 
keeping the ratio of image-lens distances: AtAB = 22h“^, magnihcation ratio 3.2.) If one 
chooses instead to shift the lens to a position that gives the observed relative magnihcation, 
the time delay is AtAB = 5.3h~^ days, but then the images would both be at a similar 
distance from the lens, which is not as observed. 

If instead we model the system as a SIS-|-constant density disk, we introduce two new 
parameters, surface density and disk length. The farther image’s light passes through the 
disk only 8.7h~^ kpc from its center, so we can start by assuming that the disk length is 
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greater than this and remove one parameter. The eqnations for a disk with radius greater 
than 8.7/i“^ become 


CTt = 


2So 

Scr 1 + 7o 

2So Oy 

Scr 7o(l + 7o) 
and the lensing potential of the disk in this case is 
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( 6 ) 


Solving for the image positions in Table I gives Sq = 7.8 x 10®/iMQ/kpc^, ay = 135 km/s, 
and a source position of (0.037, 0.12). The magnihcation ratio is 1.3, and the time delay is 
now AtAB = GSh~^ days. See Figure 1. 


The addition of the disk component has four main effects. First it breaks the spherical 
symmetry, allowing the lens not to be collinear with the images. Second, the velocity 
dispersion in the halo becomes signihcantly less, as the disk component contributes 
substantially to the deflection angle. And third, in this case the magnification ratio is 
reduced, though in general it depends strongly on the orientation of the images around the 
lens axis. Also the time delay is increased, because the potential now has a term quadratic 
in 9, and because the geometric time delays do not cancel out as they do in the SIS case. 


It is interesting to note that the surface density of the disk comes out as a reasonable 
number, although a priori it could have been anything. The deflection caused by a disk with 
density So ~ 10^h~^M q/ kpc^ is of the same scale as that caused by a halo with ay ~ 100 
km/s. Since we observe roughly this amount of luminous matter in the disk, it is clear that 
the disk’s contribution cannot be ignored. 

Of course a disk with constant surface density out to 9h~^ kpc is not a realistic model. 
We can realistically assume the disk’s surface density falls off as an exponential with 



scale length Rg. In Figure 2 the images from an exponential disk with Rg = 5h~^ kpc, 

So = 1.9 X 10®hMQ/kpc^ and a halo with = 120 km/s are shown. The images are within 
the position errors for the lens center, and the magnihcation ratio is roughly correct at 0.97. 

Interestingly we hnd that the magnihcation ratio is strongly dependent on the relative 
orientation of the images and the lens axis. Figure 3 shows the magnihcation ratio as a 
function of angle from the y-axis. A magnihcation ratio of 1.3 requires that the source 
be about 20° oh the y-axis, which also corresponds to the image orientation found in 
B1600+434. 

In this preliminary study we have found a surprisingly large disk-to-halo ratio in order 
to ht the image positions. For the exponential disk, the total disk mass is 3 x 10 ^^/i“^Mq, 
compared to a halo velocity dispersion of = 120 km/s. In comparison, standard lore 
(Binney and Tremaine 1987) gives values of 6 x 10 ^'^Mq for the disk mass and 155 km/s 
velocity dispersion for the Milky Way with a disk scale length of 3.5 kpc. While there are 
large uncertainties in these parameters for the Milky Way, our calculated disk mass does 
seem excessively large. 

Without an understanding of the uncertainties, it is hard to judge how signihcant 
these results are. In particular the inclination we have used of 65° (based on the axis ratio 
measured by Jaunsen and Hjorth 1997) could be much higher, which tends to decrease the 
scale length and hence the total disk mass. A preliminary inspection of HST archival images 
of B1600-I-434 reveals a dust lane nearly bifurcating the disk, suggesting an inclination 
greater than 80°. However, if the large disk mass did turn out to be signihcant, it may be 
evidence that the dark matter halo is battened (see for example Oiling 1996; Sackett 1996; 
Rix 1995), and the non-spherically situated dark matter is contributing to what we have 
assumed is a visible matter disk. It is also possible that the dark matter halo may have a 
core radius, and since the lensing is mostly sensitive to the matter inside the images, the 
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we are measuring is lower than the value one would use to characterize the galaxy. We are 
preparing to address these issues in a more detailed study to follow. 

4. Other effects of Including a Disk Component 

The inclination of the disk can have a dramatic effect on the lens cross section for 
multiple lensing of background sources, as we show in Figure 4. Here we model the disk 
with constant surface density again, but we choose parameters more representative of an 
L*-type spiral. 

In this case it is possible to obtain a simple expression for the critical inclination, 

7 c = cos9c, at which the cross section begins to exceed that of the SIS halo alone. Just 
note that this occurs when the inner caustic due to the disk touches the SIS caustic. If 
the inclination slightly exceeds 9c, then we can have a source just outside the SIS caustic 
and just inside the disk caustic (see the 75°-curve in Figure 4). A source at this location, 
{Px,Py) = (tts/s, 0) with asis = (As/-Ds) 47 r((Tt,/c)^, has three images. Using the lens 
equation 

(5 = 9 — ddisk — dsis (7) 

we see that, assuming the disk radius is larger than asis, one image is at 

9^o = ‘2asis/ (1 - ^ —r) and 6 ^^ = 0 . ( 8 ) 

/ \ Ucc(l + 7c)/ 

Because we know that as the source passes through the caustic the images will merge, we 
know that the other two images must be close to the hrst one. Thus there are two images 
at 6 * 3 , Ri 6 * 3,0 and ± 6 *yo, where 6 * 3,0 3> 9yo- For these two images the lens equation requires 

o^sis _ 2 So 

9x0 ^crbc(l T be) 


1 


(9) 
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independently of 6yo. Therefore, 



gives the critical inclination. 

In Figure 4 we show two inclinations for a disk of mass M^isk = and radius 

Rdisk = kpc, inside a SIS halo with = 155 km/s. For this choice of parameters 
6c ~ 73°. Thus, at an inclination of 75° the cross section shows a modest increase. By 
contrast, at 85° we hnd a very large increase by a factor of ~ 3.5 over that of the SIS halo 
alone. The disk model is not realistic, of course, but the effect is quantitatively similar for 
an exponential disk. In particular, dc is still fairly well approximated by Eq. (10). 

Whether this increase in cross section signihcantly increases the fraction of spiral 
galaxies responsible for multiple image lensing depends on the amount of mass in the 
disk. In Figure 5 one sees that for parameters like those quoted for the Milky Way the 
average (integrated over cos6i) increase is small, roughly ~ 25%. While, for a heavy disk 
like our model of B1600+434, the average cross section increases by ~ 90%. But since 
the amplihcation also increases with disk inclination, the brightest quasars lensed by 
disk galaxies will probably have highly inclined disks, especially if the disks are massive. 
Thus the fraction of lenses that are spirals, and their relative inclinations, will place 
strong constraints on the total disk mass. We shall discuss these issue in more detail in a 
subsequent publication. 

The last problem we study here is how our results would affect the gravitational 
lensing effects on DBAS. Bartelmann and Loeb (1996) and Smette (1995) have pointed out 
that gravitational lensing could signihcantly affect the observed properties of DBAS for 
two reasons. First, magnihcation bias signihcantly enhances the probability of observing 
damped absorption in a QSO spectrum. Second, there is a cutoh in the observed column 



density due to light bending, which makes it far more likely to observe absorption that 
occurs at large radii. Since the column density of absorbing HI falls off exponentially, this 
signihcantly affects the observed column densities. 

In their studies, Smette et al. (1997), Bartelmann and Loeb (1996) treated the spiral 
disk as essentially massless for lensing purposes. We hnd that including the disk lensing will 
signihcantly enhance both these effects, as we show in Figure 6. There we show the effects 
of adding disk lensing on the magnihcation and impact parameters of two sources. First, 
note that there is a signihcant increase in the magnihcation of the image that, precisely, 
occurs on the surface of the disk. For the example shown, the L*-type disk considered 
above, the point-source magnihcation increase is ~ 5.2(3.0) for the source at position 1 
(2). Second, there is an increase in the impact parameter, which would result in further 
reduction of the observed column density. These results are typical of what we hnd at high 
inclinations of 60°- 80°. Note that the increased magnihcation by inclined disks increases 
the likelihood of lensing in this conhguration; this could obviate the need for rather thick 
disks (e.g., Wolfe 1997) to account for the metal line systems seen in DLAS. Also one must 
be aware that unresolved multiple images may pass through the same disk at various places 
complicating the reconstruction of the absorbing system. 

A complete analysis is now far more involved, as one must model the disk as an 
exponential disk, average over the observed magnihcation, and average over the disk 
inclination. Also the ehects of dust in obscuring images (as clearly demonstrated in the 
relative reddening of image B to A in B1600-I-434) must be included. We will discuss the 
full problem in a subsequent publication. 
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5. Conclusions 


We have shown that the inclusion of a disk component in a spiral gravitational lens 
model has signihcant effects that cannot be ignored. Our results imply that a host of issues 
ranging from the ability of spirals to lens QSO’s to detailed properties of lens systems 
should be reconsidered to properly include the lensing effects of their disks. We have shown, 
as an example, that a SIS+disk model can account for the observations of the B1600+434 
system with fairly reasonable parameters. The implications of our results are that a velocity 
dispersion signihcantly smaller than the SIS prediction of 200 km/s, and a time delay 
greater than one month, are to be expected. 


With more data and more detailed analysis, the modeling of multiply-imaged systems 
lensed by spiral galaxies will constrain the relative amounts of matter in the disk and halo 
when disk lensing is included. However, because of the dependence on new parameters, 
further study will be necessary before such systems can give an unambiguous measurement 
of the Hubble parameter h. We have also noted that the disk contribution to lensing will 
make the amplihcation and by-pass effects on the observed distribution of DLAS in bright 
QSO’s even more pronounced at high column densities. 


After the hrst version of this paper was posted on the preprint archive as |astro-| 
ph/9701lT0| , the preprint |astro-ph/9702078l was posted by Wang and Turner discussing 
critical curves and caustics of inclined, constant density disks. Their equations and hgures 
can be seen to be the same as ours if one is careful to scale appropriately by the Hubble 
constant h, and to notice that they have calculated angular diameter distances with all the 
matter in the universe in clumps while we have assumed it to be distributed smoothly. 


This work has been supported by NSF grant PHY-9600239 at UM-St. Louis, and 
by NSF grant PHY-9402455 and a NASA ATP grant at UCSC. AM also acknowledges 
GAANN support. The authors are grateful for helpful suggestions from the referee. 
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Fig. 1.— Lensing model for B1600+434 by a SIS + uniform disk, with disk length > 8.7h“^ 
kpc. The observed images are the black dots. The small contours show the images of a small 
circular source located at the star symbol in the source plane. Since lensing conserves surface 
brightness the relative sizes of the images is indicative of their magnihcation ratio. The model 
magnihcation ratio calculated from point sources is 1.3. The source position is chosen such 
that the model images are centered on the observed images. The model magnihcation ratio 
is roughly equal to 1.3. Here = 135 km/s and Sq = 7.8 x 10®/iMQ/kpc^. The dashed 
curve is the SIS “caustic”. The model lens has two caustics: a “tangential” caustic (dotted 
line), and a “radial” caustic that coincides with the SIS caustic. Also the critical curve 
corresponding to the “tangential” caustic is shown as the outer solid curve. 

Fig. 2.— Lensing model for B1600+434 by SIS + exponential disk. The observed images 
are the black dots. The images are only slightly off their observed positions and have roughly 
the correct magnihcation ratio 0.97. Rg = 5/i“^kpc, So = 1.9 x 10®/iMQ/kpc^ and = 135 
km/s. Caustics are as in Figure 1. 

Fig. 3.— Dependence of magnihcation on source orientation. The magnihcation ratio is 
plotted vs. angle from the y-axis. The source is at a constant distance of 07126. The lens 
parameters are the same as in Figure 1. 

Fig. 4.— Dependence of multiple-lensing cross section on disk inclination. Here the lens has 
a halo with = 155h km/s and a uniform density disk with mass Mdisk = IO^^/i^^Mq and 
radius Rdisk = kpc. The dashed circle represents the caustic of the SIS halo. When 
inclined at 75° the “tangential” caustic is the thick line and the “radial” caustic coincides 
with the dashed circle. Inclined at 85° the “tangential” caustic is the thin line and the 
“radial” caustic is the dotted line that joins the dashed circle. (A source behind the tip of 
the arrow-like region would have 6 images.) 



16 


Fig. 5.— Increase in multiple image cross section to SIS alone as a function of inclination 
for a SIS + exponential disk model. The solid line is for a massive disk like that found for 
B1600+434, see Figure 2. Broken line is for a less massive disk like a Milky Way type galaxy 
with Mdisk = 6 X and radius Rdisk = S.5h~^ kpc. 

Fig. 6.— Dependence of magnihcation and bending angle on disk inclination. The solid 
lines are image contours of two circular sources (1 and 2) at positions marked by the star 
symbols. Here the disk is as in Figure 4, inclined at 75°. The dashed contours mark where 
the image contours he for the SIS halo alone. The dotted line marks the edge of the constant 
density disk. See discussion in the text. 
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Table 1. Image and Lens Positions 


Object 

Ox 

9y 

Lens 

O'/O 

0"0 

A 

-0"08 

1"12 

B 

0"26 

-0''23 


















